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5 ■ ABSTRACT 

(N 
o , 

1^ . We present results of deep direct imaging of the radial velocity (RV) planet- 

(-■h ; host star 14 Her (=GJ 614, HD 145675), obtained in the L' band with the Clio-2 

^ ! camera and the MMT adaptive optics system. This star has one confirmed 

planet and an unconfirmed outer companion, suggested by residuals in the RV 
P_( ! data. The orbital parameters of the unconfirmed object are not well constrained 

since many mass/semimajor axis configurations can fit the available data. The 
! star has been directly imaged several times, but none of the campaigns has ruled 

^ , out sub-stellar companions. With ~ 2.5 hrs of integration, we r ule out at 5cr 

?-j ! confid ence > 18 Mj companions beyond ~ 25 AU, based on the iBaraffe et al. 

^ ! (120031) COND mass-luminosity models. Combining our detection limits with 



fits to the RV data and analytic dynamical analysis, we constrain the orbital 
parameters of 14 Her c to be: 3 < m/Mj < 42, 7 < a/AU < 25, and e < 
0.5. A wealth of information can be obtained from RV/direct imaging overlap, 
^ \ especially with deep imaging as this work shows. The collaboration between RV 

\ and direct imaging will become more important in the coming years as the phase 

Csl ! space probed by each technique converges. Future studies involving RV/imaging 

^ ' overlap should be sure to consider the effects of a potential planet's projected 

separation, as quoting limits assuming face-on orientation will be misleading. 

^ ! Subject headings: techniques: high angular resolution — stars: individual (14 

! Her) — planetary systems 
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Introduction 



In recent years, several exoplanets 


lave b 


2008; 


Kalas et al. 


2008; 


Laeranee et al. 


2010) 



by far discovered the most exoplanets (525 as of February 2011, http://exoplanet.eu), the 
few direct imaging discoveries have expanded our knowledge of exoplanets significantly. The 



directly-imaged planet orbiting Fomalhaut ( iKalas et al.l 120081 ) shows evidence for having 
cleared out dusty material in the star's debris d i sk. The four ~ 1 Jupiter- mass (Mj) 



planets in the HR 8799 system (IMarois et al.ll2008t iMarois et al.ll2010l ) orbit at such a wide 
range of distances that they challenge all current planet formation theories. The sun-like star 
G J 758 has a ~ 40 Mj companion with a Tf.f f < 600 K, mak ing it one of the coldest objects 



ever discovered ( iThalmann et al.l (120091); ICurrie et a. 



helped probe exoplanet at niospheres; 



(120 101) 1. Recently direct imaging has 



Janson et al.l ( 120101 ) characterized the atmosphere of 



HR 8799 c via spectra and iHinz et al.l ( I2010l ) characterized the atmospheres of HR 8799 b, 
c, and d via photometry at 3-5 ixm. For t he first time significant orbital motion has be en 
imaged for the planet in the /3 Pic system (ILagrange et al.l (l2010l ). iLagrange et al.l ( 120091 ) ). 
This has helped constrain the planet's orbit. 

While the RV technique is mostly sensitive to planets orbiting close (semimajor axis 
a < 5 AU) to their host stars, direct imaging probes the outer (a ^ 5 AU) regions. Thus 
direct imaging and RV complement each other. A star's planetary system architecture can 
be fully determined when it is studied by both RV and direct imaging. A planet's true mass, 
semimajor axis, eccentricity, and orbitial inclination can be determined when it is detected 
by both RV and direc t imaging. Mass-luminosity models for low -mass brown dwarfs and 



gas-giant planets (e.g, iBaraffe et al.l ( 120031 ); iBurrows et al.l (120031 )) can also be constrained 
and improved in this case. 

The optimal targets for direct imaging are young planets on wide orbits. Since the 
optimal RV targets are old, quiescent stars with close-in planets, the sample of systems that 
satisfy both techniques' requirements is currently very small. Even if one cannot currently 
image any known companions, perhaps im aging of RV stars can detect previously undiscov- 
ered companions. iRodigas fc Hinzl (120091 ) used Monte Carlo simulations of RV data to show 
that as many as 15% of RV systems that contain a single, moderately-eccentric planet may 
have an additional massive planet on a wide orbit. In this case the outer planet's RV signal 
is weak enough that it is not detected above the noise. A more favorable case would be if 
the outer planet was suggested by a long-term trend in the RV data. Then we know the 
companion is there, and it is just a matter of detecting it. 

The star 14 Her (=GJ 614, H P 145675) is a prirn e target for RV/direct imaging overlap. 
At a distance of 17.6 pc ± 0.1 pc (Ivan Leeuwenll2007l ). it is close enough that direct imaging 
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can probe as close as ~ 9 AU from the star. 14 Her is a multiple-body system. It has one 
detected planet (of minimum mass 4.64 Mj, semimajor axis a = 2.77 AU and eccentricity 
e = 0.369) and a sec ond unconfirmed companion that has been s ugges t ed by a long-term 



trend in the RV data (IWittenmyer et al.l (l2007| ) jGozdziewski et al.l ( 120061 ) jGozdziewski et al. 



( I2OO8I )). Two-planet Ke plerian fits to the data s uggest that the outer c ompanion is > 2.1 
Mj an d orbits at > 7 AU (IWittenmyer et al. I l2007h . Dynamical analysis by lGozdziewski et al.l 
(120081 ) suggests a best-fit, minimum mass of ~ 8 Mj and a semimajor axis of ~ 9 AU 



though there were a host of lo w-y^ solutions. Fi g . [T] sh ows the current RV data excluding 
the orbit of the primary from IWittenmyer et al.l (120071 ). Due to the lack of data covering 
the orbit of the outer companion, there are many different possible mass/semimajor axis 
solutions. 

The companion, at a minimum, orbits on a fairly wide orbit. This is beneficial for direct 
imaging since luminous objects on wide orbits are easier to detect via imaging than objects 
on small orbits. This comes from imaging typically having better contrast (and therefore 
sensitivity) farther from the star. Imaging, especially at the current epoch, is also favorable 
given that for most of the possible orbital solutions 14 Her c is likely to have a large projected 
separation; only for some high inclination solutions can the planet be behind or in front of 
the star. 

Direct imaging of stars with known exoplan ets has yielde d stron g constraints on plan - 



etary masses orb i ting a t large separations (e.g., iJanson et al.l ( 120091 ): iHeinze et al.l ( 120081 ): 
Kenworthy et al.l ( 120091 )). With regard to 14 Her, there have been several direct imaging 
campaigns. Two direct imaging campaigns with the Keck and Lick Observatories have al- 
ready ruled out stellar-mass (M > 80 Mj) companion s beyond 9 and 12.7 A U, respectively 
( iLuhman fc Jayawardhanal l2002l : iPatience et al.ll2002l ).^ ICarson et al.l ( l2009l ) used the Palo 
mar t elescope to rule out ~ 70 Mj companions beyond ~ 18 AU.^ Recently iLeconte et al. 



(I2OIOI ) used the Advanced Electro- Optical System telescope to conduct an imaging survey 
of nearby solar-like stars, one of which was 14 Her. However with only ~ 20 minutes of 
integration they do not report any significant detection limits for this star. Deeper imaging 
is required to set meaningful constraints on 14 Her c and any other potential companions. 

To investigate the nature of 14 Her c and probe for additional high-mass planets and 
brown dwarfs, we have carried out deep direct imaging of 14 Her in the L' band with the 
MMT adaptive optics (AO) system. In Section 2 we describe the observations and data 
reduction. In Section 3 we present our contraints on 14 Her c's mass, combining analysis of 
the published RV data with our direct imaging results and analytic dynamical analysis. In 



''These limits assume a face-on orientation to the star, which is unhkely. 
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Section 4 we summarize and conclude. 



Observations and Data Reduction 



Observations were carried out at the MMT on Mount Hopkins in Arizona on the night 



of UT 2010 May 30. We used aio-2 (IFreed et al.ll2004J : ISivanandam et aLll2006l : iHinz et aL 



20061 : ICurrie et al.ll2010[ ). Arizona's high-contrast near-infrared camera, and observed in the 
L' band. The field of view was approximately 9" by 30", with a plate scale of 29.9 mas/pixel, 
determined by observations of the binary star HD 223718. We turned the instrument rotator 
off so that the field of view would rotate thr oughout the observat ions; this is essential for high- 
contrast angular differential imaging (API. iMarois et al.l (120061 )). Observing conditions were 
optimal, with clear skies, good seeing, and the AO providing consistently good atmospheric 
correction. Throughout the observations, we nodded the telescope every few minutes by 
10" so that each image would contain the target star and sky background. We obtained 
8764.7 seconds (2 hrs 26 mins) of integration on 14 Her. At ~ 50% efficiency with Clio-2 
in the L' band, this translated into ~ 5 hours of wall-clock time. A small fraction of the 
obtained images were not used in the data reduction due to the AO being off, wind shake 
distorting the star, or an infrequent occurrence of losing the star on separate telescope nods. 
We also obtained un saturated images of 14 Her and HD 203856 (an L' standard star from 
Leggett et al.l (|2003[ )) using a neutral density filter to calibrate our photometry. 



Images were saved and written as stacked data cubes. All data reduction was performed 
with custom scripts in Matlab. We divided each image by the number of frames in each data 
cube and by the integration time to produce units of dn/s for each pixel. We corrected the 
images for bad pixels using a bad pixel mask, and we removed detector and image artifacts 
as follows: since each image serves as both a target image and a background image, each 
image was subtracted from the opposite- nod image obtained closest in time to it, as long as 
that image had not already been used in the nod-subtraction process. Once an image is used 
to remove sky background, it is no longer available for nod-subtraction. Since we saturated 
the central 0."3 core of the star in each image, we determined the pixel location of the star 
by smoothing each image with a 25 pixel (0."75) disk (pillbox) average, finding the maxi- 
mum pixel location, then calculating the center of light at this pixel within a 0."75 radius. 
Calculating the center of light allows for sub-pixel registration, which increases contrast. We 
used the center of light location to register each nod-subtracted im age to a common pixe l 
location. We then reduced the images using the LOCI algorithm ( iLafreniere et al.l 120071 ). 
Each image was rotated clockwise by 2.53° minus its parallactic angle to obtain north up 
and east left. The 2.53° rotational constant was determined by observing the binary star 
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HD 223718. The final image was produced by median-combining the set of all the images. 
The change in parallactic angle between the first and last images taken during the night was 
158°, allowing us to detect point sources all the way up to the saturated central star (0."3). 



3. Results and Discussion 

Fig. [2] shows our final reduced image of 14 Her. No candidate companions are identified. 
Some point-sources masquerade as "real," but after dividing up the data into the first and 
second half of the night, all sources can be ruled out as speckles. 

To determine our sensitivity level, we calculated the standard deviation per pixel in a 
5 pixel (0."15 = FWHM) annulus centered on the star, from 0."3 to 2". We compared the 
total flux in a 5-pixel aperture centered on the unsaturated image of HD 203856 to the 14 
Her standard deviation per pixel x the square root of the number of pixels in the 5-pixel 
aperture. We used this ratio to calculate the apparent magnitude of the 14 Her background 
as a function of separation from the star. We used unsaturated images of 14 Her to calculate 
its apparent magnitude relative to the standard star HD 203856. We calculated this value to 
be 4.76 mag. We used this value and the 14 Her background annulus calculation to determine 
the contrast as a function of separation from 14 Her. Based on these calculations, we achieve 
excellent imaging contrast and sensitivity with our observations, reaching contrasts of ~ 10 
A mags at 0."4 and > 13 A mags beyond 1". 



3.1. Phase Space Constraints 



The true mass of a directly imaged luminous object depends heavily on its age. The 
younger the host star, the brighter the object is expected to be. Knowing the host star's age 
accurately is essential for pinning down the companion mass range. Unfortunately 14 Her's 
age, like many stars' ages, is fairly uncertain. Kinematically the star is a metal-rich thin- 
disk star, sugge sting an age < 10 Gyr. T h e chro mospheric activity and age-rotation-activity 
estimates from iMamajek fc Hillenbrandl (120081 ) estimate an age of 7.8 Gyr and > 8 Gyr, 
respectively. G iven that 14 Her is very metal r ich ([Fe/H] ~ 0.3), age-rotation-activity cali- 



brations (as in lMamajek fc Hillenbrandl (|2008[ )) could o yerestimate the age , since these rely 



on solar-metallicity calibrator stars. On the other hand, iRocha-Pinto et al.l ( 120001 ) estimated 
an age of ~ 3.3 Gyr from chromospheric activity and metallicity studies. However this value 
may underestimate the true age since it is based on a linear age-activity fit, which does no t 
appear to accurately fit samples with well-determined ages (IMamajek &: Hillenbrandl l2008l ). 
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There are in total s even reported age values of 14 Her. One of these is almost twice the 
age of the universe (ITakeda et al.l 120071 ). so we ignore this value. Because our observational 
result is a non-detection, we want to quote conservative upper hmits to mass and should thus 
use a conservative age in calculations. Therefore we take the maximum of the remaining 
age values^ and adopt this value, 8 Gyr, as the age of 14 Her. The COND age tracks do 
not fall on 8 Gyr, but rather between 5 and 10 Gyr. Therefore we quote Mj sensitivities 
interpolated between these two ages. 



As a first step towa r ds con s training 14 H e r c's 



phase space, we fit the available RV data 



from IWittenmver et al.l fl2007f ). iButler et al.l ([20061), and iNaef et al.l t004 \ at fixed outer 
planet mass and semimajor axis, allowing all other system parameters to vary. In Fig. [3] 
we plot the Ax^ = 1 contours for inclination angles ranging from 11° to 90° (colored lines). 
14 Her, like mos t stars with RV pla nets, has an unknown i nclina tion angle. Nonetheless 
Han et all J200lh found i > 25° and iReffert fc QuirrenbachI J201lh found 11° < i < 154°, 
both using Hipparcos astrometry. The latter range is equivalent to 11° < i < 90° since we 
are dealing with the amplitude of sini. To be conservative we adopt this range of inclination 
angles when computing 14 Her c phase space. 

A second constraint co r nes fr om analytic dynamical analysis. We used equations derived 
by ISzebehely fc McKenzid (Il98ll ) to calculate the maximum mass 14 Her c could have at a 
given semimajor axis such that the three-body system is stable. This is shown as the grey 
curve in Fig. [31 The equations assume that both planets have e = and treat 14 Her b 
as a massless test particle. We know that this is not true for 14 Her b; it is both fairly 
massive (> 4.64 Mj) and eccentric (e = 0.369), both of which would limit the outer planet's 
minimum semimajor axis. This means that our results from this analysis are conservative 
estimates. Any phase space to the left of this dynamical curve is considered unstable. 

In addition to ruling out mass and semimajor axis values, the dynamical curve also helps 
constrain eccentricity. We computed eccentricity contours from the RV fits, but do not plot 
them (for clarity). The e > 0.5 solutions are ruled out since the regions of these eccentricity 
contours that lie within the Ax^ = 1 contours correspond to dynamically unstable massive 
companions. Therefore we take e = 0.5 as an upper limit on 14 Her c's eccentricity. 

Constraints from our direct imaging cannot simply be plotted in Fig. [3] because this 
would assume a face-on orientation. For imaging we must deal with projected separation 
instead of semimajor axis. Therefore the most probable orbital elements, which are within 
the regions to the right of the dynamical curve and within the RV contours, must be mapped 



^3.33 Gyr 



Rocha-Pinto et al.ll2000l). 5.0 Gvr dValenti fc FischalboOSh. 5.24 Gvr dTakedalbOQTh. 6.9 Gyr 



20041 ). 7.8 Gyr |Mamaiek fc Hillenbrandlboosi) . > 8 Gvr |Maniaiek fc HillenbrandH20oi ) 



( Wright et al 
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into projected separation space at the epoch of the imaging observations. 

In Fig. Hlwe plot the dynamically-bound RV contours, mapped into projected separa- 
tion, for 11° < z < 90° (colored contours). At high inclinations, some orbital solutions place 
14 Her c behind or in front of the star due to projection effects. The final constraint comes 
from our direct imaging sensitivity (solid black curve in Fig. H]). The planet's allowed mass 
and separation values must reside below the imaging curve and inside the dynamically-bound 
RV contours. We shade this region grey. 

14 Her c then has the following constraint on mass: 3 < m/Mj < 42. Using 42 Mj as 
an upper limit on mass, we can constrain 14 Her c's semimajor axis. This is shown in Fig. 
ISl which is the same as Fig. [3] except that we have included this mass constraint (horizontal 
black line). The phase space above the horizontal line is ruled out. The allowed phase space, 
which is below the mass constraint line and to the right of the dynamical curve, is shaded 
grey. 14 Her c's semima jor axis is then given by 7 < a/AU < 25. Our limits agree well with 



the dynamical results of iGozdziewski et al.l (120081 ) . who estimated a best-fit mass, semimajor 



axis, and eccentricity of ~ 8 Mj, ~ 9 AU, and ~ 0.2, respectively. 



4. Conclusions: What is 14 Her c? 



A main goal in exoplanet imaging studies is to determine how planets form. Without 
better constraints on mass, semimajor axis, and eccentricity we are unable to comment on 
whether 14 Her c formed by core accretion or disk instability. However, the close separation 
and low mass constraints presented here indicate the object formed out of the disk, rather 
than via cloud fragmentation. If 14 Her c's semimajor axis is closer to ~ 7 AU, its mass is 
likely to be small and comparable to that of 14 Her b, which would suggest formation by core 
accretion. It may take many more years of additional RV and imaging observations before 
more powerful constraints can shed light on how the two planets formed. Nonetheless this 
work is important for showing the potential of RV/direct imaging overlap. We have demon- 
strated the "proof of concept," showing how much orbital phase space can be constrained 
with just a few hours of observations combined with RV data and simple dynamical anal- 
ysis. We also showed the importance of considering projected separation of an RV planet 
when constraining phase space. Future cases of RV/imaging overlap should work in pro- 
jected separation, as assuming face-on orientations can be misleading. The next generation 
of large ground- and space-based telescopes will be able to probe lower-mass regimes and will 
therefore set more powerful constraints. This deeper imaging, combined with continued RV 
monitoring of systems like 14 Her, will help characterize the architectures of the planetary 
systems beyond our own. 
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Fig. 1 . — RV data excluding the primary (b) component for 14 Her, from IWittenmyer et al. 
(120071 ). Because the observation baseline is shorter than the companion's long period, there 
are many possible good fits to the data and therefore many possible values for 14 Her c's 
mass, semimajor axis, and eccentricity. 
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arcseconds 

Fig. 2. — Final reduced image of 14 Her, obtained using the LOCI algorithm. North is up 
and east is to the left. The central 0."32 has been masked out since this region is saturated. 
The stretch is -5(T to 5(7, where a was computed as the standard deviation in 5 pixel-wide 
(0."15) annuli centered on the star. No candidate companions are identified in the image. 
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semimajor axis (AU) 



Fig. 3. — Mass vs. semimajor axis for 14 Her c, from Ax^ = 1 RV contours (colored 
lines, 11° < i < 90°) and an analytic dynamical co nstraint (grey line ) . The RV contours 



were c ompu ted by fitting the available RV data from IWittenmyer et al.l (120071 1 . iButler et al. 



( 120061 ). and iNaef et al.l ( I2004J ) at fixed outer planet mass and semimajor axis, allowing all 
other syst em parameters to vary. The dynamical constraint, which comes from equations 
derived by ISzebehely fc McKenzid ( ll98ll ) , represents the maximum mass 14 Her c could have 
at a given semimajor axis such that the three-body system is stable. 
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projected Eseparatioii (AU) 



Fig. 4. — 14 Her c's mass vs. projected separation at the epoch of our imaging observations. 
The sohd black curve is our imaging sensitivity curve. Any phase space above this hne is 
ruled out at 5a confidence. The colored contours represent the projected dynamically-bound 
Ax^ = 1 RV contours. At high inclinations, some orbital solutions place 14 Her c behind or 
in front of the star due to projection effects. The region shaded grey represents the planet's 
most probable mass and separation values. 
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Fig. 5. — The same as Fig. |3l except the 42 Mj upper hmit on mass (horizontal black line), 
which was calculated from Fig. HI is included. The phase space above this line is ruled out. 
The allowed phase space, which is below this line and to the right of the dynamical curve, 
is shaded grey. 14 Her c's maximum semimajor axis is then ~ 25 AU. 



